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The behavior of a driven granular gas in a container consisting of M connected compartments 
is studied employing a microscopic kinetic model. After obtaining the governing equations for the 
occupation numbers and the granular temperatures of each compartment we consider the various 
dynamical regimes. The system displays interesting analogies with the ordering processes of phase 
separating mixtures quenched below the their critical point. In particular, we show that below a 
certain value of the driving intensity the populations of the various compartments become unequal 
and the system forms clusters. Such a phenomenon is not instantaneous, but is characterized 
by a time scale, r, which follows a Vogel-Vulcher exponential behavior. On the other hand, the 
reverse phenomenon which involves the "evaporation" of a cluster due to the driving force is also 
characterized by a second time scale which diverges at the limit of stability of the cluster. 
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I. INTRODUCTION 



One of the most ubiquitous phenomena in nature is the existence of matter in different phases. Water, for instance, 
can form either a sohd, a liquid or a gas. When such a system is brought, by a sudden change of a control parameter, 
from a one phase equilibrium state , to a point in the phase diagram where it can exist in two different phases, it will 
reach its new thermodynamic equilibrium only through a non equilibrium process, called phase separation 1]. In the 
present paper, we stress that granular fluids possess a phenomenology which recalls closely that of standard matter. 
More specifically, we show that some processes observed during an ordinary condensation-evaporation process have a 
counterpart in vibrated granular fluids. In particular, the formation of dense patches surrounded by rarefied regions 
is similar to the phase separation dynamics associated with a liquid-gas transition 0. 

The idea that granular fluids |3|-0> i-^- large collections of inelastic particles fluidized by the action of an external 
driving force, under appropriate conditions may exhibit behaviors typical of ordinary fluids is contained in some 
recent papers. Sunthar and Kumaran reported the coexistence of different states in vibro-fluidized granular beds 
[^. Argentina et al., instead, claimed that a vibrated gas of inelastic particles displays Van der Waals loops Q in a 
pressure-density diagram. 

Some years ago, Schlicting and Nordmeier in a seminal paper considered an assembly of steel balls in a vertical 
container of height L partitioned in two connected sections, by a dividing wall of height I < L. They observed that, 
when the container was vigorously shaken, the number of balls in the two sections was statistically identical, whereas 
the two populations were dramatically different for weak shaking. The reason for such a behavior is the competition 
between the particle diffusion induced by the shaking and the tendency to cluster resulting from the inelastic collisions. 
Various authors contributed with theoretical explanations of such a problem. These range from phcnomenological 
flux models ^SJ-,I2J, to urn models .12J-J^], to kinetic theory approaches [l^. . Moreover, the Twente group 

|l2l | presented new sets of experiments which have stimulated more interest in the problem. 

We shall consider a granular gas subjected to a vigorous shaking and initially equi-partitioned into several identical 
compartments and show that it presents a phenomenology resembling that of spinodal decomposition 20]. At some 
instant t = the shaking intensity is decreased and the system evolves towards a new statistically steady state. The 
system has two possibilities: one is to persist in the homogeneous state, the second to cluster. The crossover between 
the two regimes occurs at a particular value of the driving intensity through the amplification of long-wavelength 
fluctuations. However, the analogy with spinodal decomposition is not complete, since in the late stage the order 
parameter does not saturate and the behavior of the granular system shows substantial differences with respect to 
familiar ordering systems. 

The remainder of the paper is organized as follows. In section II we specify the model, which is based on the 
assumption that the grains are described as inelastic hard disks subjected to a stochastic driving force. It is then 
possible to derive within the framework of the Boltzmann equation the governing law for the occupation number 
and the granular temperature of each compartment, the relevant variables in the problem. The governing equations 
represent an extension of those recently employed in ref. and validated by Monte Carlo simulation. In section III 
we perform the stability analysis of the homogeneous asymptotic solutions of the governing equations and extract 
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predictions about the crossover from the homogeneous regime to the clustered regime. In section IV we illustrate the 
numerical results for one and two dimensional systems. Finally, in section V we present our conclusions. 



II. MODEL 

We propose a simple extension of the model employed elsewhere in order to study the steady state properties 
of a vibro-fluidized granular gas to the case where the total volume available to the particles is divided into a series of 
identical compartments and the grains can move from one to the other by jumping over a vertical wall. An assembly of 
Af inelastic hard-spheres moves in a d-dimensional domain partitioned into M identical regions of volume, V separated 
by vertical obstacles. Each compartment contains Ni particles so that X^fii — 

When two particles collide their velocities after the collision, denoted with a prime, are obtained in terms of the 
(unprimed) precoUisional velocities through the relation 

v'l = Vl - i(l + a)(vi - V2) -6-)^ (1) 

where & is the unit vector directed from particle 1 to particle 2, and a is the coefficient of restitution. A particle in 
the i-th box, besides colliding inelastically with remaining {Ni — 1) within the same box is subjected to the action of 
a white noise random force, which compensates the energy losses due to dissipative forces. 

The dynamics of the k-th particle between two successive collisions is based on a Langevin type equation of motion 
for each grain with a fluctuating random force accounting for the action of the external driving force: 

-rr = + ^k, (2) 

at n 

where —r^^Xk is a viscous term and a Gaussian random acceleration, whose average is zero and variance satisfies 
a fluctuation-dissipation relation: 

{^kf^mnAt')) = 2— 4™<5p.5(i - t'), (3) 
rriTb 

where is proportional to the intensity of the driving [2ll | and /i, v denote vector components. The rate at which 
the kinetic energy is dissipated by collisions is proportional to 1 — a^. 

Finally, the particles contained in compartment i can migrate into compartments i ± 1 with a probability per unit 
time, T~^, provided their kinetic energy exceeds the fixed threshold Tg, which is related to the gravitational energy 
necessary to overcome the vertical barrier. 

Instead of considering the individual trajectories of the particles, one can study the single particle phase-space 
distribution function /(r, v, t), which contains most of the relevant statistical information about the system. In order 
to obtain /(r, v, t) one still has to solve the associated Boltzmann equation for /(r, v, t) 0|, which is not a simple task. 
However, a simpler description can be achieved by assuming that the relevant properties of the system are described 
in terms of the average particle population and of the average kinetic energy in each compartment. Therefore, we 
introduce the following coarse grained distributions: 



/^(v,t) = :^^dr/(r,v. 



where the integration domain is restricted to the volume of the z-th compartment. Such an approximation clearly 
neglects the gradients which are present within each box. 

Interestingly, also the coarse grained distributions /i(v,i) can be obtained by means of a Boltzmann like kinetic 
approach. In fact, fi(y,t) evolves in time due to: a) the interaction with the heat bath b) the collisions between the 
particles belonging to the same box, c) particle diffusion from one compartment to the other. The effect due to a) is 
represented by a Fokker-Planck term: 

1 ^(^A+vi)/.(vi,t) (4) 



Tb l9vi TO 9vi 

whereas the effect of b) is encapsulated in the collision term, I{fi, fi) (see ref. Finally, we add a term taking 

into account the flux of particles between neighboring compartments. Since only fast particles contribute to the latter 
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process such a term contains a Heaviside 6 function. Collecting terms we obtain the following equation of motion for 

T /' 5 \ ^ 1 ^ 1 
aj,(vi,i) =/(/„/,) + — 7^ /.(vi,<) + -T^(vi/,(vi,i)) 0(|vi|-w,)[2/,(vi,t)-/,+i(vi,<)-/,_i(vi,i)] 

(5) 

A more detailed treatment of such transport equation can be found in reference [IJl together with • 

In order to characterize the macrostate of the system we only need the average number, Ni(t), of particles in 
compartment i at instant t, and the corresponding granular temperature Ti. These quantities are related to fi by the 
equations 

(6) 



(7) 



N,{t) = J dv/,(v,i) 
where and Vi is the volume of the same compartment, and 



N, 



Wd /dvmvV.(v,t). 



The problem can be reduced to a simple set of governing equations for the occupation numbers and the temperatures 
of the various compartment (see ref. [Tg for details) if one assumes a Gaussian shape for the distribution function fi 
and d — 2 



(8) 



In the case of non communicating compartments (tj — > oo), each containing N* = J\f /M particles, the granular 
temperatures Ti assume the value T* , determined by the solution of the equation 



l + na{l - a^) 



N* T 



T, 



2V V m 

In the general case, after simple manipulations, we obtain the following set of coupled equations 



(9) 



dt r. 



dt Ts 

+ (iV,:+ie-^=/^'+i + 7V,_ie-^=/^-i - 2N,e-^^/^^)i2Ts - T,)] - 2ju;,N,T, + -7V,(Tb - T,) 

n 

where the local dissipation rate is 



(10a) 
(10b) 



2V 



(11) 



Let us notice that the ansatz (jSJ is only dictated by technical convenience and not necessary, since it leads to a small 
set of coupled equations for Ni and Ti. In fact, one way to improve systematically the approximation is to assume 
fi to be given by a Gaussian multiplied by a linear combination of orthogonal polynomials, the so called Sonine 
polynomials. However, in this case one should consider extra variables besides Ni and Ti and the resulting set of 
equations would be of higher rank. 

We turn now to consider the most relevant behaviors associated with the governing equations (|10a|l and Ijl0b|l . 
Incidentally we notice that eqs. (|10a|l and l)10b|l in the case of AI = 2 are identical to the corresponding equations of 
reference [l2f after substituting Ts with 2ts- 

A large collection of grains initially distributed uniformly in an array of M identical compartments may remain 
homogeneous or not according to the intensity of the driving at which is subjected. The Twente group observed 
experimentally a homogeneous configuration under vigorous shaking, and clusters for weak shaking The crossover 
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from the first to the second regime occurs when the "thermally" induced diffusion is not sufficient to prevent the 
spontaneous tendency of the grains to form clusters, due to the collisional cooling. Such an instability is the counterpart 
of the separation process which occurs in a system consisting of two compartments only. In such a case the left-right 
symmetry, i.e. the difference between the left population, Ni^, and the right population, iV/j, is spontaneously broken 
below a certain critical temperature. Fig. la illustrates the corresponding behavior of the asymmetry parameter 
A = \Nl — Nji\/J\f versus the total number of particles for a fixed value of Ti, — 0.7, whereas in Fig. lb we display the 
variation of the granular temperatures within the two compartments versus the total number of particles. The curve 
bifurcates at a critical values of the total number of particles identifying a temperature, Tc, below which the left-right 
symmetry is broken. We shall discuss how such a mechanism manifests in the case of many compartments and gives 
rise to a phenomenology similar to that of the spinodal decomposition In general the space of dimensionless 

control parameter is large since the system properties are functions of a, M, J\f, V/a'^, Th/Tg, Tb/Ts and ^Ts- In the 
following we shall measure temperatures in units of Tg, areas in units of and time in units of Tb/2. 




2000 4000 6000 8000 10000 2000 4000 6000 8000 10000 

PARTICLE NUMBER PARTICLE NUMBER 

FIG. 1: Phase diagram relative to a system of two compartments: in Fig. la we show the asymmetry parameter A = 
\Nl — Nr\/N as a function of the total number of particles, TV. In Fig. lb granular temperatures of the two compartments 
are displayed as a function of N ■ The heat-bath temperature is fixed at Th = 0.7, Ts = 1, rj, = 2 and Ts = 0.5, while the 
compartment volume is 100 and a — 0.7. The granular temperature is measured in non dimensional units. 



III. STABILITY ANALYSIS OF THE HOMOGENEOUS STATE 



In the case of M identical compartments with cyclic boundary conditions, the choice Ni = N* and Ti = T*, where 
T* and N* are related by eq. Q represents a uniform solution of eqs. (|10|l . for all values of the control parameters. 
We observe that T* is the granular temperature of a system of M particles, equally distributed into M compartments 
of volume V and subjected to a heat bath Ti,. 

On the other hand, it turns out that such a uniform solution is stable only at high temperature, where a diffusive 
mechanism tends to restore any small perturbation about the homogeneous state. On the contrary, the uniform state 
below a certain temperature turns out to be unstable with respect to spontaneous fluctuations, due to the clustering 
mechanism induced by inelasticity. As shown in Fig. lb, the granular temperature in the case of a system with 
two compartments takes on two different values when the total population exceeds a threshold value. In the case of 
many connected compartments a related phenomenon occurs. In order to illustrate it, we introduce a small amplitude 
sinusoidal perturbation T/ = T* + 5Tk ex-p{ikl), and Ni ~ N* + SNk exp{ikl), where k — 2TTn/M, with n = 1, .., Af — 1 
and / = 1, ..,M denotes the compartment. 

Expanding linearly eqs. H10a() and IjlObj) about the symmetric fixed point T*,N* one finds the result: 



SNk = ~- 2(1 - cos(fc)) 



N*T, 



2 



(r*)2 

] (1 - cos(fc)) + ( 37W* 



STk 



-{T* + 2Ts){l - cos(fc)) + juj*T* 



(12a) 



(12b) 
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FIG. 2: Variation of the eigenvalue A„(fe) with respect to the wave- vector k for different choices of the heat-bath temperature 
Tj, = 0.28 (continuous hne), Ti, = 0.38 (dashed line) and Tt — 0.45 (dotted line). The data refer to a system consisting of 
100 compartments, each initially filled with 180 particles. The compartment volume is 100 and a — 0.7. The threshold 
temperature is Ts = 1, Tt = 2, Ts — 0.5: Notice that the mode becomes unstable, i.e. X„{k) becomes positive, below the critical 
temperature. 



The associated eigenvalues A„(fc) and Xxik) of the dynamical matrix of coefficients correspond to the two relaxation 
modes of the system. The larger eigenvalue A„(fc) vanishes quadratically when A; — > 0, reflecting the conservation of 
the global number of particles. However, due to the coupling between the density and the thermal fluctuations, for 
finite k, A„(fc) presents a variety of behaviors, associated with different physical phenomena. 

In order to classify these behaviors, we consider the following expansion valid for small k values 



(13) 



Above Tc, a2 is negative, the eigenvalue A„(fc) describes a diffusion process by which a local density fluctuation is 
re-adsorbed. In other words, the particle fluxes caused by the coupling to the external driving are sufficient to restore 
homogeneity. 

Below a certain characteristic temperature, Tc, a fluctuation which increases locally the population is amplified. 



The local granular temperature drops due to the increased collision rate, since (from eq. 0) ui oc NiT^ (x , and 
Ti(xNi ' . Thus the particles arriving from the other compartments remain trapped, causing a further reduction of 
the local temperature. This phenomenon is described by the formula 



a2 



Tc 

— - 1 



i.e., a2 becomes positive below Tc which is given by 



(14) 



(15) 



with juj* ^ (t(1 — a^)|^^-^. Contrary to the previous diffusive case, particles tend to cluster. 

According to the sign of 04 (for a2 > 0) the initial regime is different. In fact, if 04 < , X{k) may display a 
maximum at a finite wave- vector, km < 27r (T;, = 0.38, in fig. El , whereas for a4 > 0, A„(fc) attains its maximum only 
in correspondence of the largest wave-vector for Ti, — 0.28. 

Correspondingly, in the second case the growth process initially resembles the early stage of the spinodal decompo- 
sition process, during which the homogeneous state is unstable with respect to long wavelength fluctuations. These 
are exponentially amplified, whereas short scale fluctuations decay. 

Physically speaking, a local density increment induces a decrement in the granular temperature. If the temperature 
of the bath is sufficiently low or the number of particles sufficiently large this leads to an instability, i.e. more particles 
will flow toward the region where the temperature is lower. Diffusion will not be able to compensate such a tendency. 

Finally, there exist a second collective mode Xrik), which essentially describes how temperature fluctuations decay. 
It is always negative, due to the presence of dissipation caused by friction and collisions, and given by 



Xrik) — — S'-fUj 



ck^ 



(16) 
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IV. NUMERICAL RESULTS 
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FIG. 3: Evolution of the occupation profile relative to a system constituted by 100 compartments and 180 particles initially 
placed in each compartment. The temperature, is chosen below the critical temperature and is Ti, = 0.38, while the remaining 
parameters are the same as in fig.|21 The three snapshots from left to right refer to times ti = 500, t2 ~ 1500 and ts — 3000 in 
units of Tb- 
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FIG. 4: Growth of the profile width relative to a system constituted by 100 compartments and 180 particles initially placed in 
each compartment. The heat-bath temperature is Tb = 0.38 and the remaining parameters are the same as in fig.|H The right 
figure represents the parametric plot W{t) versus the average granular temperature for the same system. 
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FIG. 5: Time evolution of the average number of occupied compartments, f{t), relative to a system with M = 100 com- 
partments and 180 particles per compartment at t — 0. The curves refer from right to left to granular temperatures 
Tg = 0.273,0.26,0.25,0.24,0.23 (Fig. 5a). The remaining parameters are those of fig El In Fig. 5b we display the varia- 
tion of the average granular temperature for the same system (Fig. 5b) . 



In the present section we shall study numerically the properties of the model. We have integrated the eqs. H1U|) 
with a = 0.7, V = lOOcr^, = I, Th — 2, = 0.5. We have chosen a as unit length, as energy unit and 2ts as 
time unit. The value of Tb and of TV and AI varied from case to case, together with the initial conditions. 
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FIG. 6; Characteristic time associated with the position of the kink, shown in Fig. 5a, versus average granular temperature. 
The system is a linear array with M = 100 compartments. The remaining control parameters are the same as in fig|5| The 
behavior is not Ahrrenius-like. The data can be fitted by a Vogel-Fulcher law r„/ = ^exp[A/(Tg — To)], where the fitting 
parameters are A = 200, A = 0.072 and To = 0.293. 

The stability analysis discussed so far, describes only the linear, early stage of the evolution. In order to explore the 
behavior in the late non-linear regime we have solved numerically the governing equations Hl()|) . We considered a one 
dimensional array of compartments, initially equally populated, at the same granular temperature, i.e. Ni — N* and 
Ti = T* and added a small random perturbation. According to the previous linear analysis, two different behaviors 
can be observed as the control parameters, such as the heat bath temperature and the average density, vary. At fixed 
N, for large values of Tb, such that Tg > Tc, the initial perturbation is re-adsorbed diffusively, while at small values 
of Tf, {Tg < Tc) the perturbation is exponentially amplified. In the latter case, the coUisional cooling determines 
a decrease of the local temperature in correspondence of the regions more populated and clustering begins. Some 
compartments, randomly selected by the dynamics, act as germs for the nucleation process illustrated in figs. 01 
After the initial regime few compartments grow at the expense of the remaining which become empty. The distance 
between highly populated compartments increases, since they compete for particles. However, unlike the late stage 
spinodal decomposition process ji^l, we observe that the domains do not grow in size, but in height. This feature can 
be understood because no saturation mechanism is present in the model, so that the occupation number in a single 
compartment can become of the order of the total population Af. The growth process occurs by diffusion against the 
density gradient, since particles move from low populated regions toward highly populated regions. The evolution of 
the width, W{t) (see fig. EJ 



illustrates quantitatively how the process occurs. During the initial stage W{t) remains much smaller than N* and 
only after a characteristic time r it begins to rise steeply. In fig. ^ we show W{t) versus time and the parametric plot 
W{t) against the average value of the granular temperature. We notice that the latter plot recalls an order parameter 
versus temperature plot in a system undergoing a phase transition. 

A second quantitative measure of the clustering phenomenon is represented by the following statistical indicator 



The "entropy" h is non negative, vanishes when all particles are confined in a single compartment and takes on its 
maximum value, ln(M), when all compartments are identically populated. Thus / — exjp{h) represents a measure of 
the number of occupied compartments. In fig. [51 we display the evolution of the average granular temperature and 
of / in various situations. Above Tc the indicator / relaxes toward M, whereas in the low temperature region, due 
to clustering, / decreases toward a plateau value P < M. Interestingly, such a relaxation time t increases as the 
system approaches the temperature Tc from below. The temperature dependence of r close to T^ displayed in fig. 
is consistent with the Vogel-Fulcher law 




(17) 




(18) 



T,f=Aexp[A/{T-To)]. 



(19) 



The dependence of the characteristic time t„/ on the temperature is a direct consequence of eqs. (|13|l and (|14|) . We 
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also remark that the plateau value of /, reached by the system for t > Tyf increases as Ti, decreases, indicating that 
the system remains trapped in some metastable configurations. 

We turn, now, attention to a different process obtained by considering the evolution of an initial configuration, in 
which all the particles are located inside a single compartment at t = 0. According to the level of one can observe 
two different processes: a) for large Tj, the occupation number in the central cell decays toward the fully symmetric 
state Ni = Af/M and limt_>oo f{t) — M; b) for small the occupation of the compartment remains constant. 

In fig. [71 we display the variation in time of the population in the compartment, where it was initially placed, for 
various values of the heat bath temperature. We observe that N{t) decreases more and more slowly as the transition 
temperature is approached from above. 
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FIG. 7: Variation of the height of a cluster, initially containing 1200 particles, as a function of time for temperatures Tj, — 0.30, 
0.35, 0.40, 0.45, 0.50, 0.55 and 0.60 from top to bottom. 

In fig. IHlis shown the temperature dependence of the average time r necessary to wash-out the initial single cluster 
configuration. We observe that r = C/{Ti, — Tp)"^/^ diverges at the crossover temperature Tp. In particular the plot 
of the occupation number versus time indicates the growth of a plateau when Th Tp The smaller the temperature 
deviation from the limit of stability Tp the longer the plateau. 

In fig. Owe display the transition temperature at which a cluster of = 1200 particles placed in a cell at i = 
"explodes" as a function of the number of boxes, M. Such an effect has been observed experimentally by Lohse et al. 
|l0l | and explained in terms of their flux model. 
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FIG. 8: Characteristic time associated with the explosion of a cluster of 1200 particles located at t = in a single compartment 
versus the heath bath temperature. The system consists of a linear array of 100 compartments, while the remaining parameters 
are the same as fig. |5| 

Finally, we consider the inverse process, by which a cluster instead of "evaporating" , grows at the expenses of its 
neighbors. Such a study also provides some information about the interfacial properties of the model. We prepare 
the system in an initial configuration where the population Nc of a single compartment is much larger than that of 
the remaining compartments. If the profile were uniform with Ni — Noo — 180 and — 0.7 the state would be 
linearly stable with respect to perturbations. However, if we place Nc = 1800 particles in one of the compartments 
and keep the remaining populations at Ni = 180 we observe the following scenario: the spike grows in height, while 
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NUMBER OF COMPARTMENTS 

FIG. 9; Stability of an isolated cluster as a function of the number of compartments. We show the phase boundary between 
the stable and the unstable regime relative to a cluster of A'^ = 1200 particles initially placed in a single compartment. The 
vertical axis represents the heat-bath temperature at which the transition takes place, while the horizontal axis represents the 
number of compartments. 



the compartments closer to it slowly empty. As shown in fig. 1101 the population profile develops a gradient due to 
the flux of particles from the "bulk" to the "tower" . Its shape is similar to the depletion layer associated with the 
growth of a liquid droplet in a "sea" of over-saturated gas [23|. In addition, we observed that since the flux terms 
in eq. H12a() are very small the local values of the granular temperature and of the occupation number are related 
by eq. jnj. The profile varies almost linearly indicating that the process can be assimilated to a diffusion controlled 
interfacial growth. 




FIG. 10: Evolution of the depletion layer around a growing cluster at times t = 500, t = 2500 and t = 20000, for = 100. 
The initial configuration is A''^ — 180, A^so = 1800 and Tb = 0.7. The system starts with a central peak and a supersaturated 
"bulk". As the system evolves we observe the increasing thickness of the depletion layer. Matter accumulates in the central 
compartment. 



It is worthwhile to comment that the model discussed above presents a shortcoming, since the local density can 
grow unbounded until all the particles in the system occupy the same well. In fact, the approximation scheme does 
not treat adequately the mutual repulsion between the particles at high densities. Particles are allowed to pile up 
and reach arbitrarily large packing fractions. In order to eliminate such a shortcoming of the model, it is possible to 
introduce phenomcnologically a new ingredient into the theory. The model was modified introducing an occupation 
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number dependence in the characteristic time r,: 

Ts{N^)^To{l-N,/Mf'^ (20) 

That is the probabihty of escaping from a compartment becomes larger when it is occupied by a large number of 
grains. This correction is based on the results of microscopic studies [2^. The form above has been chosen in order 
to preserve the form of the equations as simple as possible and by phenomenological requirements. We observed that 
during the the early stages of the process, where the linear analysis of section III applies, there are no appreciable 
changes to the situation described in the previous sections. On the other hand, the late stage presents a different 
scenario. When the compartments saturate the occupation number and the granular temperature become locally 
stationary and the growth proceeds in the neighboring compartments and so on. 

V. CONCLUSIONS 

To summarize, we introduced a model for compartmentalized driven granular gases and studied it using the methods 
of kinetic theory. We have found a rather rich "phase" behavior and the emergence of new qualitative properties as 
the number of particles becomes sufficiently large. We have pointed out that the system undergoes a long-wave length 
instability and orders in a fashion similar to the process which occurs during the spinodal decomposition in fluid 
mixtures. However, the late stage of the process is radically different, because the granular gas does not possess a 
surface tension mechanism which restores homogeneity. Thus the usual competition between bulk and surface free 
energy cost which determines the growth of larger and larger domains is not at work here. 

The present approach, in spite of its simplicity, has the advantage with respect to the so called flux models of 
relating the microscopic parameters to the macroscopic observables in a natural fashion. It can be extended to treat 
granular mixtures, where the granular temperature of each component must be treated as an independent variable 
to be determined self-consistently 113, [l^ ■ I'^ addition, the approach can be improved by including non Gaussian 
corrections to the distribution function or by solving numerically by the Direct Monte Carlo simulation method the 
Boltzmann equation. 
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